Goos-Hanchen-Like Shifts in Atom Optics 
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We consider the propagation of a matter wavepacket of two-level atoms through a square potential 
created by a super-Gaussian laser beam. We explore the matter wave analog of Goos-Hanchen shift 
within the framework of atom optics where the roles of atom and light is exchanged with respect 
to conventional optics. Using a vector theory, where atoms are treated as particles possessing two 
internal spin components, we show that not only large negative but also large positive Goos-Hanchen 
shifts can occur in the reflected atomic beam. 
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I. INTRODUCTION 



In conventional optics for light waves, Goos-Hanchen 
in 1947 discovered that a light beam under the condition 
of total reflection can experience a lateral shift (or dis- 
placement) along the surface of a dielectric boundary 
This pioneering-work has stimulated a large volume of 

studies Jl UMM EMM [jlMJMI^^ 

concerningthe Goos-Hanchen shift of reflected or trans- 
mitted 0,13) 111 light beams of different polarizations [HI 
in different media characterized with, for example, peri- 
odic structures [3, [H [H, [3(1 , (left or right) handedness 
0, HI 11, |3, multilayers 0, weakly absorbing JTp, HI] , 
lower-index [1 91 ] or negative- index of refraction [la ] . etc. 
The key physics behind the Goos-Hanchen shift is the na- 
ture of wave interference. From the perspective of wave 
optics, the incident beam of a finite transverse width can 
be viewed as composed of plane wave components, each 
of which has a slightly different transverse wavevector. 
Each wave component, after the total internal reflection, 
undergoes a different phase shift, and the superposition 
of all the reflected wave components gives rise to the lat- 
eral shift of the intensity peak in the reflected beam Q ■ 



In this sense, it is not so much the total internal re- 
flection but rather the phase modulations for different 
plane wave components that remains the true mecha- 
nism behind the lateral shift. Thus, the Goos-Hanchen 
shift is expected to occur in matter waves where particles 
have finite masses. As is known, under the usual condi- 
tions (or temperatures), electrons possess a de Broglie's 
wavelength much longer than atoms because the latter is 
much heavier in mass than the former. Thus, it is much 
easier to demonstrate the Goos-Hanchen shift with elec- 
trons [H, HII or even neutrons [H, |36j] than with atoms. 
The situation, however, has been rapidly changed over 
the last two decades. Nowadays, ultracold atoms with a 
relatively long de Broglie's wavelength can be routinely 
made available, thanks to the rapid advancement of the 
laser cooling and trapping technology. Motivated by 
the fact that ultracold atoms have led to many impor- 
tant applications in atom optics [37}, we explore, in this 



paper, the matter wave analog of Goos-Hanchen effect 
within the framework of atom optics where matter waves 
of ultracold atoms are manipulated by laser fields. An 
important difference between the matter waves in atom 
optics and the light waves in conventional optics is that 
atoms have internal electronic structures while photons 
are structureless. Thus, an accurate description of the 
Goos-Hanchen effect in atom optics must regard atoms as 
particles possessing internal spins (energy states). To the 
best of our knowledge, our work here represents the first 
that is seriously devoted to the problem of Goos-Hanchen 
effect with cold atoms. As such, we limit our goals to 
establishing a general theoretical framework and to ap- 
plying it for a basic understanding of the Goos-Hanchen 
effect in matter waves with cold two-level atoms, while at 
the same time hoping that our work can draw significant 
attentions from experimentalists for future applications. 

Our paper is organized as follows. In Sec. [Til we derive 
a set of coupled 1-D Schrodinger equations to describe the 
scattering of two-level atoms by a super-Gaussian laser 
beam in a 2-D setting. In the same section, we present 
the connection between the lateral shifts and the coeffi- 
cients of reflection and transmission. In Sec. IIII1 we 
derive, with the help of a dressed state picture, a set of 
analytical expressions for the reflection and transmission 
coefficients, which are to be used in Sec. IIVI to signifi- 
cantly simplify our calculations. In Sec. IIVI we combine 
the tools developed in Sec. [H] with those in Sec. IIIII to 
numerically investigate, within the context of atom op- 
tics, the matter wave analog of Goos-Hanchen-like shifts. 
Finally, a conclusion will be given in Sec. [V] 



II. MODEL AND BASIC EQUATIONS 

Figure[T]is the schematic of our model, in which a mat- 
ter wavepacket composed of two-level atoms of transition 
frequency tu a is obliquely incident upon a "slab" made up 
of a travelling laser beam of frequency u>l and wavenum- 
ber kh- In our model, we require that both the atom and 
laser beams be sufficiently wide along the direction nor- 
mal to the plane of incidence (x — y plane) so that the 
degree of freedom in the z-dimension can be completely 
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FIG. 1: A wave packet of atoms with two internal states 
impinges on a laser "slab" . 



liberated. Under such a circumstance, we can adopt the 
following coupled 2 — D Schrodinger equations [38| 



• nv 2 T n(x) _ ikrVT 

2m 2 
W 2 



«*2 = - 



-A-i- * 2 - 
2m 2 ' - 



(la) 
*l> (lb) 



to describe the evolution of the wavefunctions, *f?i and 
\I/2, of the ground state |1) and the excited state |2). In 
Eqs. |T|), we have defined m as the atomic mass, V 2 = 
d 2 /dx 2 + d 2 /dy 2 as the 2-D Laplacian operator, A = 
wl — uj a as the laser detuning, and 7 as the decay rate 
of the excited atomic state. In addition, we describe the 
dipole interaction between the laser field and atoms by a 
Rabi frequency in the form of 



Q(x) = QF(x) , 



(2) 



where f2 is the peak value and F(x) is a normalized 
spatial function, representing the laser profile. In this 
paper, the laser is assumed to possess a beam profile 
in the form of a high-order Gaussian function, F (x) = 



exp 



(x 2 /w\ 



,N 



; such a beam is experimentally acces- 
sible through spatial shaping techniques [39|, H(J El[ or 
optical techniques [42}. To further simplify the prob- 
lem, we restrict our study to the super-Gaussian beam 
with an order number N so large that, to a fairly good 
approximation, F (x) can be idealized as a step function 

F(x) = ( 1 ' -V2<*<172, (3) 
w [0, x> L/2,x< -L/2, y ' 

Next, we utilize the fact of the Rabi frequency [Eq. 
P])] being y independent to eliminate y in favor of y- 
wavevector k y through the Fourier transformation 



#1 (r, t) = J dkyfa (x, k y , t) e Lk yy-^\ (4a) 
* 2 (r, t) = [ dk y fo {x, k y , t) eWv+^y-^*, (4b) 



where r = {x,y). By doing so, we transform Eqs. ([T]) 
into coupled 1-D Schrodinger equations 



h 2 d 2 

2m dx 2 



l + V)4> 



(5) 



where <j) — {(pi, ^2) is a two-component vector field, / is 
a 2 x 2 unit matrix, and V is the potential matrix given 
by 



V = — 



n 
n 2(<5 + i 



n ) 



In Eq. 



we have defined 



Hk y k L 



2m ' 



(6) 



(7) 



as the effective detuning, where HkykL/rn and hk 2 L /2m 
are, respectively, the Doppler and the photon recoil fre- 
quency. Equation ([5]) serves as the starting point for next 
section, where we calculate the scattering matrix, which 
determines all the scattering properties of our model. 

For now, we turn our attention to the lateral shifts 
of the reflected and transmitted wavepackets. For this 
purpose, let's first jump ahead to Eqs. (fT5j) of Sec. 
IIII1 which define various transmission and reflection co- 
efficients via the stationary scattering solutions in free 
space. Of relevance to our interest here are the coeffi- 
cients of transmission T\ and reflection R\ of the ground 
state; we ignore T2 and R2 of the excited state since 
the excited wave, being highly susceptible to the decay 
by the spontaneous emission, cannot propagate far from 
the scattering region. Let 0f ,T (k) be the phases of re- 
flection and transmission coefficients defined through the 
relation 



Si (k) = I Si (k) 



^ef(k) 



(8) 



where for notational simplicity, we have used (and will 
continue to use) S = R and S = T to symbolize reflection 
and transmission, respectively. For an incident wave 
packet initially (t = 0) located at (xn, yo) far away from 
the laser beam (see Fig. [T]), we can construct, through 
the superposition of the time-independent solutions [Eqs. 
(|15p]. its reflected and transmitted wavepackets in the 
form of 



*f (r, t) = J d 2 k'f (k' - k) I Si (k') |e J "i S ( k '' r ) (9) 
where ■q l : arc the total phases defined as 

(k',r) = Gf (k') - k' x Ax R + k' y Ay R - (10a) 

ril (k',r) = ej 00 + k' x Ax T + k' y Ay T - Of*, (10b) 

" 2m 



3 



with Ax R = x + xq, Ax t = x — xq, Ay R — Ay T = y — yo, 
and k' 2 = k' x 2 + k' y 2 . In Eqs. ©, /(k'-k) is a 
(real) weighting function peaked around k' = k with a 
momentum distribution sufficiently narrow and smooth 
that Eqs. ([9]) represent fairly accurately the plane waves 
of velocity Kk/m. Equation ((9]) indicates that the re- 
flected and transmitted waves are the result of interfer- 
ence among different wave components distinguished by 
wavevector k' = (k' x ,k'y). As a result, the values of 
these waves at a given time and location, (r,t), depend 
crucially on the phases [Eqs. (fTT)|) ] of each k' component. 
In particular, fy RT (r, t) reach peak values when a con- 
structive interference takes place or equivalently when 
Vk'?7i ' at k = k vanish. Using this condition, we 
find that the peaks of the wavepackets in the coordinate 
space propagate with time according to 



Ax R 
Ax T 



.— £t+— L, Ay R 
m ok r 



hk 



t ■ 



dej 

dk„ 



m 

hk 



m 



dk y 

del 

dk„ ' 



(11a) 
(lib) 



where all the derivatives are assumed to be taken at k' = 
k. 

Let to, tfi and tr be, respectively, the time duration 
of the atomic beam between t = and right before it 
hits the boundary at x = —L/2, between t = and 
immediately after it is reflected from the boundary at x = 
—L/2, and between t = and right after it is transmitted 
from the boundary at x — L/2. In terms of to, we have 
xq = —(L/2 + Hkxto/m) and yo = —hkyto/m, which, 
when substituted into Eqs. (fTTjl for Ax R ' T , allows us to 
find t R ^ T 



Ate 



ts - to 



dGf 



(12) 



Finally, by incorporating these results into Eqs. (TTT]) 
for Ay R ' T , we find that the induced Goos-Hanchen-like 
lateral shift yn due to the reflection and yr due to the 
transmission (see Fig. [T]) are governed by 



hk v . def 
Vs = ~^At s - -^i. 
m ok v 



(13) 



where dQf/dki (i = x,y and S = R,T) are evaluated 
using 



~dk~ 



1 dSi 1 91511 



Si dh | Si | dh 



(14) 



which is a direct consequence of Eq. ([S]). In contrast 
to the usual shifts, which are solely determined by the 
part directly proportional to At$, the lateral shifts in 
Eq. (fT3]) contain an additional term def /dk y . This is a 
unique aspect of atom optics, where momentum conser- 
vation during the photon emission and absorption makes 
the effective laser detuning S [Eq. |(7J)] k y dependent. 



As a result, the phase 0f becomes a function of k y via 
its dependence on S, which, in turn, leads to a finite 

aef /dk y . 

From this derivation, it is clear that (a) the Goos- 
Hanchen-like lateral shifts are the wave phenomena, that 
depend crucially on the ability of the optical potential to 
modify the phases of various matter wave components, 
and (b) the key to the lateral shifts is the transmission 
and reflection coefficients, which will be the focus of our 
study in the next section. 



III. TRANSMISSION AND REFLECTION 
COEFFICIENTS 

In this section, we construct the reflection and trans- 
mission coefficients, starting from the stationary scatter- 
ing solutions of Eq. (|5|) for an incident ground atomic 
beam having an energy E x = Hk x /2m and wavenumbcr 
k x along the x dimension. Let's first introduce the re- 
flection and transmission coefficients for the ground state, 
i?i and Ti, and those for the excited state, i?2 and T 2 , via 
the scattering solutions in regions I and III. By virtue of 
the decoupling between the excited and ground states in 
free propagation regions I and III outside the laser slab, 
the scattering solutions take the form 



e lkix +R l e~ %kix 



R 2 e~ 
T ie lklX 

J ik 2 X 



where we have defined the free-space wavevectors 



771 / 'T 

h = k x , k 2 = \l2—[5 + i- 



kl. 



(15a) 
(15b) 

(16) 



The excited-state and ground-state components in re- 
gion II are, however, mixed because Eq. (O is a coupled 
equation. To solve Eq. and thus, to find the vector 
wavefunction (p 11 in region II, we first seek to diagonalize 
the matrix V [Eq. ©] by looking for the eigenvectors of 
V. This leads to two eigenvalues V+ and V_, given by 



(17) 



The corresponding eigenvectors are expressed as the first 
and second column vectors of the following transforma- 
tion matrix 



U = 



SXD.(f 



COS if 



-e^ cos ip e lt3 smtp 



where ip and /3, defined as 

m 



tan ip 



2 7+ 



,V+ = \V + \e i ?, 



(18) 



(19) 



4 



are two angles introduced to characterize the dressed 
states 



|+) = sin<p|l) - e l/3 cos<yj|2) , 
|-) = costal) + e - lP smip\2) . 



(20a) 
(20b) 



The inverse of [/, which will also be an important part 
of the scattering problem involving vector matter waves, 
is given by 



jj-i _ , / e 1/3 simp —cosip 
' ' e tl3 cos tp sin ip 



(21) 



where / = {eT 1 ^ sin 2 ip + e i/3 cos 2 tp) is a normalization 
factor. In the absence of spontaneous decay 7 = 0, we 
have (3 = 0, and U becomes unitary. The dressed states 
can be simplified into the ones well-known in quantum 
optics [43| ■ The presence of the spontaneous decay ren- 
ders the Hamiltonian nonhermitian [3], which is why U 
is no longer a unitary matrix when [3 ^ as one can 
easily verify from Eq. (jT8|l . 

With these preparations, we now express ifi 11 , in terms 
of wavefunctions <f>± on the dressed-state basis, as 



= U 



1 A 2 e a2X + B 2 e~ a2X 1 



(22) 



where A4 and Bi are the superposition coefficients, and 
for easy organization, we introduce a% = a+ and a 2 = 
a_, where 



a± = y / 2mV±/h 2 - fc2. 



(23) 



To facilitate the derivation bellow, besides U and U~ l , 
we also introduce the matrix W and its inverse W~ x , 
where 

or equivalently, 

Wa = (ihy 1 Ua, (W- 1 )^ = ikj (t/- 1 ) .. . (25) 

Next, we require that ip 1 ' 11 ' 111 and their derivatives be 
continuous at location x = L/2, leading to 



- < 1 T2elk2 L/2 ) « 



. x / r ie * ^ 2 

V2 J { T 2 e^ L ' 2 



(26a) 
(26b) 



where (it,, Uj) are a set of new variables defined as 

Ui = A,e a ^ +B ie - a ^, (27a) 
m = a t (A,e a ^ - B t e~ a ^ (27b) 



Similarly, application of the continuation conditions at 
location x — —L/2 results in 



»g-r )- (28a) 

where (xi,yi) are defined in terms of (Ai,Bi) as 



x % = A ie - a >- + B t e^~, 
Vi = OL i (A ie - ai i -B^ 



(29) 



By combining all these conditions (for details see the Ap- 
pendix), we arrive at a set of compact formulas for the 
transmission coefficients 



Ti = 



To 



2M. 



(+) 



2MiP 



-i(k 1 +k 2 )L/2 



and for the reflection coefficients 

1 M^M^-M^M^ ' 



R 2 



where 



MiVM^-M^MiV 



m 2 2 IV1 21 „-i(k 1 +k 2 )L/2 

& 1 



(30a) 
(30b) 

(31a) 
(31b) 



M 



(±) 



n=l,2 

cosh (a n L) — i 



nj 



sinh (a n L) 



(32) 



In the next section, Eqs. (|30a|) and (|31a|) will be used in 
connection with the results from Sec. |TT] to numerically 
determine the lateral shifts. 



IV. 



DISCUSSION 



In this section, we carry out a numerical study of the 
lateral shifts by first obtaining dOf/dki from Eq. ([5]) 
with the help of Eqs. (|30a[) and (|31ap . and then deter- 
mining the lateral shifts using Eqs. (| 1 2|) and ([13]). In 
our calculation, we replace k x = kcos9 and k y = fcsin#, 
and correspondingly, d/dk x = cosOd/dk — k^ 1 am9d/dd 
and d/dky = sindd/dk + k^ 1 cos9d/d9, where k is the 
magnitude of wavevector and 9 is the incident angle of 
the atomic wave. In addition, we adopt the following 
scaled variables: A = A/7, = fi/7, &l = fei/fe 7 , 
k = k/kry, L = L/k^ 1 , V± — V±/tvy, a± = a±/k 7 , and 
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FIG. 2: (Color online). The left column calculated under 
A = —100 contains (a) S, (b) V- , and (c) 5_ as functions of 
6. The right column calculated under A = 200 contains (d) 
S, (e) V+ , and (f) a+ as functions of 6. Other parameters 
are 7 = 1, k = 3, L = 6, = 20, and k L = 8.1125. 




FIG. 3: The left column contains (a) \Rt\ 2 , (b) 9f , and (c) 
yn as functions of 6. The right column contains (d) |Ti| 2 , 
(e) Qi , and (c) i/t as functions of 6 under the condition that 
A = —100 and the rest parameters are the same as in Fig. [2] 



Vs = Us/kj 1 , where fc 7 = ^2m^/h. In all the exam- 
ples given bellow, unless stated otherwise, 7 = 1, k = 3, 
L = 6, O = 20, and k L = 8.1125. 

Let's first consider a case where the laser detuning is 
set at A = —100. At this A, the effective laser detuning 
S remains (deeply) red detuned across all the incident 
angles as shown in Fig. [^a). Under such a circum- 
stance, we have \V+\ >> \V-\ according to Eq. (| 1 T[) and 
ip approaches a small value according to Eq. (fl"9"|) . As 
a result, the scattering properties in this case are largely 
determined by the |— ) dressed state. For this reason, 
we only display in Fig. HJb) potential V- as "seen" by 
the atoms in state |— ), which corresponds to a potential 
well since Re(V-) remains negative. As a result, the — ) 
mode function oscillates in the x dimension with a spa- 
tial frequency of Im(a_), whose value is shown in Fig. 
He)]. 

The corresponding intensity Si| , phase Of, and lat- 
eral shift ys of the reflected (5* = R, left column) and 
transmitted (S = T, right column) waves are displayed 
in Fig. [3] as functions of the incident angle 9. As our dis- 
cussion above suggests, we expect to see the phenomenon 
of resonance scattering by potential wells. Indeed, in the 
region where the phase experiences a n shift [Fig. EDJb)], 
the reflectivity |i?| 2 [Fig. OUa)] approaches a small value 
while the transmissivity \T\ 2 [Fig. [3{d)] becomes nearly 
perfect, which explains the oscillatory behavior exhibited 
both in |i?| 2 and in |T| 2 . It needs to be stressed that the 
phase of the reflected wave, as shown in Fig. [3fb), in- 



creases (or decreases) sharply, as 9 (or k x ) sweeps across 
each resonance. This causes the reflected wave to expe- 
rience a large negative Goos-Hanchen shift around each 
resonance as indicated in Fig. [3Jc). 

Next, we consider a situation where A = 200. With 
this A, 8 remains (deeply) blue detuned for all the in- 
cident angles as Fig. [2jd) illustrates. As a result, ac- 
cording to Eqs. |T7]) and HU), \V+\ « \V-\ and 99 ap- 
proaches 7r/2; the scattering behavior is dominated by 
the |+) dressed state. Here, the laser beam, as illus- 
trated in Fig. (He), creates an effective repulsive potential 
i?e(V+) > 0. Under such a circumstance and provided 
that fi 2 P/2m > Re (V+), we can introduce a critical an- 
gle defined as 9 C = cos -1 y / 2mRe (V + ) /h 2 k 2 at which 
the (x) kinetic energy of the atomic beam equals the 
height of the potential barrier. (Such a critical angle does 
not exist for the case of red detuning.) In our case here, 
we identify from Fig. [2]^f) that 9 C w 69.4°, which has the 
physical meaning that bellow 9 C , the |+) mode oscillates 
at a spatial frequency close to Im (a+) while beyond 9 C , 
it undergoes quantum tunneling with \/Re{a + ) being 
the characteristic tunneling distance. 

Indeed, Fig. [4] shows that |<Si| 2 , 6f, and ys beyond 
the critical angle 9 > 9 C are qualitatively different from 
those within the critical angle < 9 < 9 C . For 9 > 9 C , 
besides relatively sharp features around the boundaries, 
both the phase and lateral shift exhibit no oscillations. 
This description applies both to the reflected and trans- 
mitted beams. For < 9 < 9 C where the incident (x) 
kinetic energy exceeds the potential height, the reflec- 



6 



tion [Fig. IDJa)] and transmission [Fig. 0|d)] oscillate. 
The phase of the reflected wave, while still undergoes a tt 
shift, decreases (or increases) sharply, as 9 (or k x ) sweeps 
across each resonance [Fig. B^b)], a behavior completely 
opposite to the case of red detuning [Fig. G^b)]. As a 
result, we see that the reflected wave develops a large 
but positive Goos-Hanchen-like lateral shift around each 
resonance [Fig. ffljc)] (This is in contrast to the phase 
and lateral shift of the transmitted wave, which remain 
relatively monotonic within the critical angle.) 

It needs to be emphasized that large negative shifts 
around resonance have been the focus of several recent 
studies for light waves that propagate throu gh a bsorptive 
medium slabs in conventional optics [ToL Il5l l26l| . In our 
atom optics model here, the ground-state matter wave is 
coupled to the excited-state matter wave, and this cou- 
pling greatly enriches the physics concerning the lateral 
shifts. Not only do we see large negative lateral shifts as 
in Fig. [3tc) when the laser is red detuned, but also large 
positive lateral shifts as in Fig. [He) when the laser is blue 
detuned. Moreover, with our atom optics model, great 
controls over not only the position but also the linewidth 
of these peak shifts can be achieved by taking advantage 
of lasers being highly tunable both in intensity and in 
frequency (not shown). 



V. CONCLUSION 

In conclusion, we have established a theoretical frame- 
work for studying the matter wave analog of Goos- 
Hanchen-like effect in an atom optics model where a 
super-Gaussian laser beam acts as a "medium slab" for 
a matter wave of two-level atoms. We have developed 
a vector theory based upon a set of coupled Schrddinger 
equations for describing the scattering of a wavepacket of 
two-level atoms off a square potential. We have derived a 
set of analytical formulas for the transmission and reflec- 
tion coefficients, which have greatly facilitated the study 
of Goos-Hanchen effect in vector models where atoms are 
treated as particles possessing two internal spin compo- 
nents. It is important to stress that the coupling be- 
tween the ground and excited components in the vector 
model combined with the tunability offered by the laser 
field creates new opportunities for studying the lateral 
shifts. In particular, we have found that in our atom 
optics model, not only a large negative Goos-Hanchen 
shift as in conventional optics [H, [lH, H^, [33] but also a 
large positive shift can take place in the reflected atomic 
beam. 
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APPENDIX A: APPENDIX 

In this Appendix, we provide the steps leading to 
Eqs. (|30|) and (|3 1 1) . To begin with, we insert, 

Ai = \ (m + e-^/ 2 and B t = § (u* - e a > L / 2 

obtained from Eqs. (f2"7} . into Eqs. l|29p . enabling us to 
express {xi,y{) in terms of (ui, Vi) as 



Xi = Ui cosh (cxiLi) — Vi sinh (aijZ) /ai, 
Ui = Vi cosh (otiL) — onui sinh [pnV) . 



(Ala) 
(Alb) 



By combining Eq. (|28ap and Eq. (|28bp . we eliminate R\ 
and i?2 simultaneously from Eqs. (|28P and arrive at a 
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single-matrix equation where is a 2 x 2 matrix given by 



which, by virtue of Eqs. (|26p and (|A1|) . is shown to be 
equivalent to 

M<+> -C^""). (A3, 

I 



M (±) = 17 



cosh(aiL) 

cosh(a 2 L) 

cosh(aiL) 

cosh {cfrzL) 



u- 1 - U 



sinh (a±L) /a± 

sinh (oizL) /a 2 

ai sinh (aiL) / 

a2 sinh (a 2 L) 



w- 1 



(A4) 



By a similar procedure, we find from Eq. (|28|) that 

»(^)-KS)-"'(£)- (A5 » 

which, with the help of Eqs. (|26j) and (jAlj) . is equivalent 
to 

/ 2jRie ^L/2X / Tiplfcl L/2X 

^ 2R 2 e ik * L / 2 )~ M \ T 2 e^ L / 2 ) ' [Ab) 



where ' is also 2x2 matrix given by Eq. (1A4|) . A 
straightforward calculation involving the use of Eq. (1231) 

shows that the matrix element My of Eq. (|A4|) has 
a simple and explicit form given by Eq. (|32p . Finally, 
we arrive at Eqs. (150)) and (I?!]) by solving Eqs. (|A3|1 and 
(IA6|) simultaneously. 

*wpzhang@phy.ecnu. edu.cn 

f ling@rowan.edu 
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